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Abstract. Garzon et al. (1997), Lopez-Corredoira et al. (1999) and Hammersley et al. (2000) have identified in 
TMGS and DENIS data a large excess of stars at 1=27° and b=0° which might correspond to an in-plane bar. 
We compared near infrared CAIN star counts and simulations from the Besangon Model of Galaxy on 15 fields 
between 15° and 45° in longitude and -2° and 2° in latitude. Comparisons confirm the existence of an overdensity 
at longitudes lower than 27° which is inhomogeneous and decreases very steeply off the Galactic plane. The 
observed excess in the star distribution over the predicted density is even higher than 100%. Its distance from the 
sun is estimated to be lower than 6 kpc. If this overdensity corresponds to the stellar population of the bar, we 
estimate its half-length to 3.9±0.4 kpc and its angle from the Sun-center direction to 45±9 degrees. 



Key words. Galaxy: structure - Galaxy: stellar content 
1. Introduction 

Star counts have been used for years to examine at vari- 
ous levels of detail the stellar contents in the Galaxy (see 
Paul 1993), whose structural parameters of the various 
morphological components are still far from being com- 
pletely known. In the last two decades there has been 
a combined effort in this direction with the use of de- 
tailed models of stellar galactic distribution (Bahcall & 
Soneira 1980; Robin & Creze 1986; Wainscoat et al. 1992; 
Lopez-Corredoira et al. 2002) along with large area, high 
sensitivity and multicolour star count surveys. The near 
infrared (NIR) members of these surveys (Eaton et al. 
1984, Garzon et al. 1993, Hammersley et al. 1994, Lopez- 
Corredoira et al. 2000, Epchtein 1997, Skrutskie et al. 
1997) are notably useful for the analysis of the galactic 
structure because of less interstellar extinction compared 
to the optical bands, in particular in the hidden on-plane 
areas of the inner Galaxy. It is in these zones that the 
morphological structures are less well studied. 

The effect of the high extinction towards the inner 
Galaxy along the Galactic plane, added to the position 
of the Sun very near to the plane itself (Hammersley et 
al. 1995), makes it particularly difficult to detect large 
scale stellar structures in the central part of the Galaxy. 
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Galaxy: disk 

There is an increasing consensus that the Galaxy is of 
barred type (see Garzon 1999 for a review), but there is 
still a considerable controversy about the morphology and 
stellar content of the bar. 

Hammersley et al. (2000) has recently derived argu- 
ments for the Galactic bar to be of half-length of roughly 4 
kpc and with a position angle with the sun-center direction 
of around 43 deg, by observing from new NIR star count 
data an excess of stars at longitudes on the Galactic plane 
at 1<27° and assuming them to belong to the bar pop- 
ulation. To continue this research work (see also Lopez- 
Corredoira et al. 2001) in the field of the structure of the 
inner Galaxy we have made use in this paper of the large 
NIR star count database obtained with the CAIN camera 
at the Teide Observatory (Hammersley et al. 2000) to- 
gether with the Besangon model of the Galaxy (Robin & 
Creze 1986; Bienayme et al. 1987; Robin et al. 2003) to ex- 
amine the excesses in both the magnitude and the colour 
histograms of the data with respect to the model predic- 
tions. This study permits us on the one hand to confirm 
with an independent method the existence of this extra 
density, and allows us on the other hand to derive consid- 
erations about its structure. Indeed, since the Besangon 
model does not contain any galactic component other than 
the thin disc (and possibly the triaxial bulge at the inner- 
most field) in the regions of interest, as will be described 
in the next sections, any excess observed in the stellar 
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density compared to the model prediction may then be 
used to analyze the morphological extent of such an extra 
population. 

This paper is organized as follows. We briefly describe 
the star count data and the Besangon model used in 
this work. We then discuss the method of deriving the 
extinction distribution along the line of sight from the 
data alone, which is then used in the model for the pre- 
dictions. Those predictions are then compared with the 
data in two different planes: magnitude and colour counts, 
which both serve for the analysis. We provide conclusions 
describing the model itself and the derived geometry for 
the component responsible for the excesses. 

2. The Besancon Model of the Galaxy 

The approach of the Besangon model is slightly different 
from the ones used previously to study the bar. Here,we 
give a brief description of the model, dealing particularly 
with the thin disc population model. More details can be 
found in Robin & Creze (1986), Bienayme et al. (1987) 
and Robin et al. (2003) 1 . 

2.1. Global description 

The Besangon Model of Stellar Population Synthesis, de- 
veloped in the visible, near and mid infrared, aims at giv- 
ing a global 3-dimensional description of the Milky Way, 
including stellar populations such as the thin disc, outer 
bulge, thick disc and spheroid, as well as dark halo and 
interstellar matter. 

Its approach is semi-empirical: both theoretical 
schemes (stellar evolution, galactic evolution, galactic dy- 
namics) and empirical laws found in the literature are 
used. Boltzmann and Poisson equations are used to make 
the model dynamically self-consistent in the direction per- 
pendicular to the plane. Artificial catalogues of stars are 
produced, giving both the fundamental characteristics of 
stars (age, distance, absolute magnitudes, kinematics) and 
observable ones (apparent magnitudes, colours, proper 
motion, radial velocities) directly comparable with obser- 
vations. Absorption, photometric errors and Poisson noise 
are also added to make simulations as close as possible to 
observations. 

The Besangon model has been used to constrain the 
density law or the luminosity function of the stellar pop- 
ulations: thin disc (Ruphy et al. 1996; Haywood et al. 
1997), thick disc (Robin et al. 1996; Reyle & Robin 2001), 
spheroid (Robin et al. 2000) and bulge (Picaud & Robin, 
in preparation). 

2.2. Thin disc 

In the studied region, thick disc and spheroid populations 
give a negligible contribution, and those from the bulge 

1 A version of the model is available on the web at 
http:/ /www. obsbesancon.fr/fesg/modele_ang.html 
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0.15-1 


0.0268 


1-2 


0.0375 


2-3 


0.0551 


3-5 


0.0696 


5-7 


0.0785 


7-10 


0.0791 



Table 1. Axis ratios of the 6 age components of the old 
thin disc. 



likely only contaminate star counts at the innermost field 
(1=15°). Thus, the only dominant stellar population taken 
into account by the model is the thin disc. Its density is 
modeled as the Einasto law (1979): the thin disc is divided 
in 7 age components and each component distribution is 
described by an axisymmetric ellipsoid. The first compo- 
nent (0-0.15 Gyr) is called the young disc and the 6 other 
ones form the old disc, which only is presented here be- 
cause of the very small number of stars belonging to the 
young one. The density law of each ellipsoid of the old 
disc is the subtraction of two modified exponentials, the 
second one representing the central hole: 



p = po x [exp(-^0.25 + (— Y) - exp(-^0.25 + (—)*)] 

with a? = R 2 + (f ) 2 , where: 

— R and Z are the cylindric coordinates 

— e is the axis ratio of the ellipsoid which increases with 
the age of the component. The 6 axis ratios of the old 
thin disc are presented in Table l_T21 

— Rd is the global scale length of the disc and Rh is 
the global scale length of the hole. In the present ver- 
sion of the Besangon model, we use i?d=2.37 kpc and 
i?ft,=1.31 kpc, deduced from model fitting toward the 
inner Galaxy (Picaud & Robin, in preparation). 

— The normalization po was deduced from the local lumi- 
nosity function ( Jaheiss et al, private communication) . 

The distributions in M v , logT e //, logg are deduced 
from an evolution model described in Haywood et al. 
(1997). Absolute magnitudes in K are computed from 
the corresponding V absolute magnitudes using the semi- 
empirical model of atmospheres from Lejeune et al. (1997, 
1998). Disc density parameters and further changes in the 
luminosity function are explained in Robin et al. (2003). 

2.3. Extinction 

The Besangon model includes its own extinction model 
which follows the young thin disc distribution, with a scale 
length of 4 kpc, an equivalent scale height of 140 pc, and 
a normalization which can be modulated. However, this 
model can be replaced by another distribution as needed. 
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For this study, we have preferred to model the extinc- 
tion without using the Besancon simulations to avoid a 
possible bias in the comparisons of the model with the 
data. The determination of the extinction distribution in 
each field will be described in section 0] 

3. Observations 

Since 1999 the IAC group has been building a Galactic 
survey by using CAIN, the NIR camera on the 1.5-m 
Telescopio Carlos Sanchez (TCS) (Observatorio del Teide, 
Tenerife). Observations of a series of 20 x 12 arcmin 2 (i.e., 
a total area of approximately 0.07 deg 2 ) have been made 
along the Galactic plane between I = 0° and I — 220° , to- 
gether with some off-plane series (mainly in b=±2°,±5° 
and ±10°). The survey uses the J (1.25/mi), H (1.65/mi), 
and Kshort (2.17/im) standard bands, with limiting mag- 
nitudes of 17, 16.5 and 15.2 respectively. This provides 
almost one magnitude deeper coverage than 2MASS or 
DENIS surveys in the K-band (Struskie et al. 1997; 
Epchtein et al. 1997), which is very useful for the study of 
the inner galactic regions such as those considered in this 
paper. For the fields we used here, the seeing was typically 
of 1 arcsec, and data were obtained only in photometric 
conditions. The main characteristics of the 15 chosen fields 
and photometric bands used are summarized in Table 

4. Extinction from K-giants 

Extinction changes very much from one field to another 
in and close to the Galactic plane, and a good determi- 
nation of its distribution is needed to compare data and 
simulations in this region. This is why a global extinction 
model is not sufficient and the distribution of extinction 
along the line of sight must be determined field by field. 

4.1. Method of extraction of the extinction distribution 

Lopez-Corredoira et al. (2002) (hereafter L02) developed 
a method to obtain the star density and interstellar ex- 
tinction along a line of sight by extracting a well-known 
population (spectral type K2III) from the infrared colour- 
magnitude diagrams (hereafter CMD). The method is ex- 
tensively described in L02, so we give only a brief summary 
here. 

(J — K, tjik) CMDs are built for each field. In the dia- 
grams, stars of the same spectral type (which would mean 
that they have about the same absolute magnitude but 
are located at different distances from the Sun) will be 
placed at different locations in the CMD. The effect of 
distance alone shifts the stars vertically, while extinction 
by itself shifts the stars both horizontally and vertically. 
Combinations of these effects, plus some intrinsic disper- 
sion in absolute magnitude and/or spectral type, cause 
the red clump giant stars, which constitute the majority 
of the disc giants (Cohen et al. 2000; Hammersley et al. 
2000), to form a diagonal broad branch running from top 
left to bottom right on the CMDs. 




Fig. 1. CMDs for two of the fields studied here. Maxima 
obtained by a Gaussian fitting to star distribution in the 
red clump strip are over-plotted. The solid line shows the 
fitted trace assigned to the red clump population affected 
by varying extinction along the line of sight. 

In order to isolate the red clump sources, count his- 
tograms have been made by making horizontal cuts (i.e. 
running colour) through the CMDs at fixed ttik- A 
Gaussian function was then fitted to the histogram to de- 
termine the position of the peak in each cut. 

The extinction Ak (wjf), to a distance r, can be de- 
termined by tracing how the (J-K) of the peak of the red 
clump counts changes with ttik- The extinction is calcu- 
lated for any given tjik using the measured (J — K) of 
the peak, the intrinsic mean colour excess definition, and 
the interstellar extinction values for 4^ and 4^ given by 
Mathis (1990). 

(J-K) 



A K = 



(J-K)o 



1.61 



(1) 



Finally, the extinction law (Ak vs. r) can be built 
for each field by using a suitable transformation for the 
distance along the line of sight: 



10" 



: -M K +5-A K ( r ) 



(2) 



Uncertainties in the method, such the metallicity ef- 
fects for the clump giants, the effect of the dwarf contami- 
nation in the K-giant strip, or differences between the real 
red clump distribution and the Gaussian function are fully 
described in L02. 

4.2. Absolute magnitude calibration 

The method described above can provide spatial informa- 
tion from the CMDs, and the only assumption being made 
is that the absolute magnitudes of all the sources being ex- 
tracted is more or less fixed. In order to deduce the distri- 
bution of extinction along the line of sight, absolute mag- 
nitudes and colors of stars must be assigned. In L02, the 
peaks were identified as pertaining to the K2III stars since 
they are by far the most prominent population, accord- 
ing to the updated "SKY" model (Wainscoat et al. 1992; 
M. Cohen, private communication) (see Fig. 2 in L02 for 
details). The K2III population mean absolute magnitude 
and the intrinsic color were assumed to be Mk — —1.65 
and ( J — -ftf)o=0.75, with a gaussian dispersion of 0.3 mag 
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1 (deg) 


b (deg) 


area (x 10 z deg z ) 


completeness limits (band used) 


15 





6.95 


16.8 (J) 


14.8 (K s ) 


20 





6.97 


16.8 (J) 


14.8 (tf s ) 


21 





6.99 


16.8 (J) 


14.8 (if s ) 


26 





12 


16.8 (J) 


14.6 (K.) 


26 


2 


6.95 


16.4 (J) 


14.8 (Jf.) 


27 





6.97 


16.4 ( J) 


14.4 (K s ) 


27 


0.5 


6.97 


16.4 (J) 


14.8 {K s ) 


27 


-0.5 


6.96 


16.4 (J) 


14.8 (if s ) 


27 


1 


6.91 


17 (J) 


15 (K s ) 


28 





6.93 


17 (J) 


15 (if.) 


32 





6.92 


17 (J) 


15 


33 


-2 


7.03 


16.4 (J) 


14.8 (K s ) 


37 


2 


7.03 


16.8 (J) 


15.4 (if) 


40 





7.05 


16.8 (J) 


14.8 (K 3 ) 


45 





6.99 


16.6 (J) 


14.6 



Table 2. CAIN survey selected regions used in this paper. 




Fig. 2. Distribution of K absolute magnitudes and (J — 
K) colours according to the Besangon model. The K0 
counts are represented with dotted lines, the Kl with 
dotted-dashed lines and the K2 with dashed lines, while 
the continuous line represents the sum of the three counts. 
The maxima of the Kl plots correspond to the red clump, 
with an absolute magnitude of -1.85 and an intrinsic colour 
of 0.64. 



in absolute magnitude and 0.2 in colour. However, in the 
present study, the absolute magnitude and intrinsic color 
values must be the ones used in the Besancon model to 
get a consistent analysis of the data. The Besangon model 
predicts a very prominent population of Kl giants with an 
absolute magnitude of Mx=-l-85 and an intrinsic colour 
(J - A') o =0.64 (FigEJl. These values are slightly different 
from those considered initially in L02. The determination 
of spectral types in the model is taken from De Jager & 
Nieuwenhuijzen (1987). 

Differences in using either the values used by L02 or 
those of Besangon are quantified and shown in Fig. [21 They 
cannot be a critical issue in the method, since the extrac- 
tion of the peaks for the star distribution in the K-giant 
strip is independent of the subsequent analysis of the spec- 
tral type corresponding to those peaks. The change in the 




5 10 15 5 10 15 

distance (kpc) distance (kpc) 

Fig. 3. Extinction along the line of sight (Ay vs 
r) obtained for the fields 1=20°, b=0°(left) and 
1=32°, b=0° (right). The extinctions obtained assuming 
M K =-IM and (J - A) o =0.75 (L02) are shown with tri- 
angles while values with Mr-=-1.85 and (J — A)o=0.64 
(Besangon model) are shown with diamonds. Differences 
are negligible, and they are only noticeable for off-plane 
regions where the extinction is lower, hence increasing the 
relative difference. 

absolute magnitude and intrinsic colour of the population 
affects only the values extracted for Ak and r. The dif- 
ferences are less than 7% in extinction and less than 9% 
in distance (see Fig|2J). Those errors are also of the same 
order of magnitude, or even less, than the intrinsic uncer- 
tainties of the real values for the absolute magnitude and 
intrinsic colour of the red clump population. 

4.3. H band. 

The analysis is made by using (J — K^wik) CMDs for 
all the fields, except one (1=37°, b=— 2°) where the K s 
band was not available. A (J — H^mn) CMD was used 
instead. The main difference between using the H band 
or the K s band is related to the absolute magnitude and 
the intrinsic colour of the red clump population. Those 
values are not as well determined for the H filter as for 
the K s , and they present a bigger dispersion than in K s 
(Wainscoat et al. 1992). As explained in H4.2I we use the 
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values defined in the Besancon model for the red clump 
population in this band (My=-1.65, (J — H)q=0.5). The 
procedure of deriving the extinction is exactly the same as 
described before, but replacing Mk and (J — K) by Mh 
and (J — H). We are analyzing an off-plane field, and the 
extinction is low enough to produce significant differences 
in the extinction distribution obtained by using either the 
H band or the K« band. 



5. Data vs. model predictions 

This study consists of comparisons of simulations from the 
Besangon model with CAIN data at longitudes 15°-30° 
to detect and study the overdensity observed in previous 
studies (Hammersley et al. 2000, Lopez-Corredoira et al. 
1999). In this region, model counts are dominated by the 
thin disc population. This is why first comparisons be- 
tween simulations and data were needed in disc fields in 
and out of the plane to verify that the Besancon model 
reproduces well the thin disc counts. Once this verification 
has been made, counts are compared on the one hand in 
the plane varying the longitude and on the other hand at 
different latitudes to extract information about both the 
horizontal extent and the vertical thickness of the detected 
overdensity. 

Model simulations have been done assuming photo- 
metric errors close to those of the data and an extinction 
along the line of sight determined by the method described 
in section 0] 

Comparisons are made in K s magnitude and J-K s 
colour for all fields except for the field at 1=37° b=+2° 
where the K s band was not available and was replaced by 
the H band, as explained before. 

The values of limits given in Table[3]show that in fields 
with high extinction, data are not complete in the J band 
for the faintest and most reddened stars. Star counts have 
then been produced until K s =13.5 (H=14 at 1=37°) and 
J-K s =3 to ensure completeness. These new thresholds do 
not guarantee perfect completeness for all fields, but this 
study needed counts deep enough to observe the overden- 
sity. 

In order to reduce the contamination of foreground 
dwarfs in star counts, only stars with J-K s >1.5 for in- 
plane fields and J-K s >0.75 for off-plane ones 2 were se- 
lected. But the estimation (given by Fig. using the 
Besancon simulations, of the dwarf contamination shows 
that there is still up to 20% of dwarfs at faint K s and at 
low J-K s for in-plane fields, and up to 40% at low J-K s for 
off-plane ones. This contamination implies uncertainties in 
this work, especially for a quantitative study. 
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2 The reddening being low, comparisons have to be done at 
lower J-K s . This implies a substantial contamination of dwarfs 



for off-plane fields. 



Fig. 4. Percentage of dwarfs in Besangon simulations by 
bins of 3-K s (left) and K s (right), for in-plane fields (up), 
and off-plane ones (bottom) . A step between two zero lines 
(dashed lines) correspond to 50% of dwarfs. 

5.1. Disc fields 

Model counts were compared with data on 3 fields on the 
galactic plane: b=0° 1=45°, 40° and 32°, and 2 off plane 
fields: 1=37° b=+2° and 1=33° b=-2°. 

Extinction distributions, model and data colour- 
magnitude diagrams and count histograms are presented 
in Figure [5] 

Model and data CMDs and histograms are rather close 
to each other, but in some fields there is a slight ex- 
cess of reddest stars. The model of disc stellar evolution 
(Haywood et al. 1997) might overestimate the duration of 
late type stars and perhaps the oldest stars of these types 
have already become white dwarfs or planetary nebulae. A 
correction has been made by rejecting giants with a spec- 
tral type later than K5 and age greater than 7 Gyr but 
does not eliminate all the excess. The problem of late type 
stars may be enhanced in off-plane fields: in these fields, 
stars are less reddened and most of them are below the 
cuts in K,, = 13.5 and J-K s =3. Moreover, because of the 
very small reddening, histograms of these fields may be 
more contaminated by foreground dwarfs even with the 
blue cut in J-K s . 

In any case, this excess of simulated stars is small 
enough to consider that the Besancon model reproduces 
rather well the thin disc and to continue the study on 
fields at 1<30°. 
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Fig. 5. Extinction distributions along the line of sight, colour magnitude diagrams and colour histograms of inter- 
mediate disc fields. The CMD on the left and the dashed line at the histogram panels correspond to the Besancpn 
simulations, while the CMD on the right and the solid line in the histograms display the CAIN data. Star counts of 
the colour histograms (in number of stars per magnitude bin and per square degree) are taken for 11.5 < K s < 13.5 



5.2. Along the Galactic plane 

In addition to the 3 on-plane disc fields (1=32°, 1=40°, 
1=45°), we compared model and data on 6 fields at 1<28° 
and b=0° : 1=28°, 1=27°, 1=26°, 1=21°, 1=20°, 1=15°. 
Extinction distributions, colour-magnitude diagrams and 
histograms are presented in Figure |H1 

Figure presents the histograms of differences 'data 
minus model' along the galactic plane. In order to com- 
pare the relative excess at different longitudes apart from 
changes in disc density from one field to another, we com- 
puted relative differences data-model w j iere m [ s th e av _ 
erage of simulated counts over the full magnitude range. 



These histograms show clearly the existence of an ex- 
cess of stars 3 at 1<27° (except for the field at 1=21°), con- 
trary to fields at 1>28° where data and model star counts 
are similar. 

The observed excess is too large to be explained by 
a bad determination of the extinction distributions along 
the line of sight. Furthermore, color magnitude diagrams 
in the Figure [S] show that the giant clumps are rather 
well reproduced for most fields, and histograms in J-K s 
in the same Figure have similar shapes : only the heights 
of the maxima change. Moreover, data are not always to- 
tally complete for K close to 13.5, especially for the most 



3 It should be noted that the position in K 3 and J-K s of the 
overdensity depends on the extinction in each field 
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Fig. 6. Colour magnitude diagrams and colour histograms 
explanation about the content of the different panels. 



j-k, data j-k 
of fields at 1<28° and b=0°. See caption of Figure for an 



reddened stars in the fields with highest extinction, which 
may attenuate the data excess in these fields. 

A possible contamination in the counts by stars of the 
Scutum Spiral arm cannot explain the observed overden- 
sity a 1<27° as the impact of the arm is only expected at 
the position of the tangential point to the lane of stars, 
located at 1=33°, with a negligible contribution for 1<31° 
(Hammersley et al. 1994). 

Except for the discrepancies in the most reddest part of 
the CMDs (which corresponds to an excess of stars while 
the great discrepancies at low longitudes correspond to a 



deficit of the model) , comparisons of disc fields have shown 
that the Besangon model reproduces rather well the galac- 
tic disc between 1=32° and 1=45°. Thus, a bad modeling 
of disc scale length or equivalent disc scale heights cannot 
explain the excess of stars, and changes made by realistic 
modifications of them would not be sufficient to remove 
this excess. 

But this is not the same thing for the hole scale length 
which, is not relevant at 1>32°. Fig. [S] gives the results 
in 3 fields b=0° and 1=15°, 20°, 25° for 3 different hole 
scale lengths: 0.5, 1 and 1.5 kpc. Values lower than 0.5 
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Fig. 7. Difference distribution vs. magnitude along the Galactic plane. Dashed lines are zero reference lines. The left 
and center histograms present the absolute difference 'data minus model' in J-K s colour and K s magnitude, and the 
range between two zero lines corresponds to 70000 mag _1 deg -2 . The right ones present the relative difference in K s 
counts(m being the average of model counts for 11.5 < K s < 13.5) and the step between two zero lines corresponds 
to 400%. Star counts are compared in 1.5 < J-K s < 3 and 11.5 < K s < 13.5. 



kpc would not significantly change the histograms, and 
values greater that 1.5 are not very realistic. But between 
these two limits, one can see that the value of the scale 
length has a great influence on the excess. Nevertheless, 
this influence decreases with longitude and cannot explain 
all the excess. But it implies another uncertainty for a 
quantitative study. 



To summarize, the analysis of the histograms of Fig. 
confirm the existence of an overdensity with respect to 
the thin disc: predictions. The overdensity begins at 1=27° 



and extends at least until 1=15°, but seems to be inhomo- 
geneous and locally disappears at 1=21°. 

5.3. Out of the Galactic plane 

We have studied 4 fields out of the galactic plane, 
roughly at longitudes coinciding with the end of the over- 
density presented in the previous section : 1=27° b=- 
0.5°, +0. 5°, +1°, and 1=26° b=+2°. The colour-magnitude 
diagrams and the difference histograms of theses fields are 
presented respectively in Figure 151 and ITUl 
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Fig. 9. Extinction distributions, colour magnitude diagrams and colour histograms of fields at b^0°. Explanations of 
the graphs are given in Figure 



The difficulties in modeling the extinction (even if it 
is low) in these fields or to select the giant clump because 
of the low number of stars, added to the problem of the 
Bcsangon model predictions on disc fields out of the plane 
make a detailed study of the overdensity difficult out of 
the galactic plane. 

However, Figure shows that the excess decreases 
with distance from the plane. Indeed, this excess present 
at 1=27° b=0° still appears at b=-0.5°; It also appears at 
b=+0.5° but some doubt remains considering the uncer- 
tainties; at b=l° the peak is not significant; at b=2° the 
excess has disappeared. 

Thus, Figure ^| shows that the excess observed at 
b=0° on different longitudes seems to be confined very 
close to the galactic plane, at least at its end (l=26°-27°). 

5.4. Characterizing the extra density 
5.4.1. Quantitative considerations 

It is very difficult to quantify the excess density. Indeed, 
imperfections in thin disc modeling or extinction model- 
ing and completeness of data do not prevent the detection 



of the overdensity but have an effect on the determina- 
tion of its characteristics such as density and position. We 
try however to extract some information about the excess 
density from the comparisons. The most significant pa- 
rameter is the star count difference at the maximum of 
the peak in K histograms. Table l5"4~D gives the values for 
the 6 fields where the overdensity has been detected and 
the peak is rather well defined. Even if values must be 
taken with great caution, they show that the excess stars 
are at least as numerous as the disc ones. 

5.4.2. Spatial location of the overdensity 

Hammersley et al. (2000) estimated with CAIN data the 
position in the H/J-H CMD of the excess peak at 1=27° 
and deduced, using values of absolute magnitude and 
color, a distance of 5.7±0.7 kpc from the Sun for excess 
stars at this field. Using the distances for the stars from the 
Besangon model (which has other values for absolute mag- 
nitudes and colors), we have tried to replicate the same 
determination in the K s band making use of the same 
CAIN database. 
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Fig. 10. Histograms of differences. Explanations of the graphs are given in Figure Hand the same scales are used. 
However we have changed the cuts on J-K for histograms due to the fact that the reddening is now lower than in the 
in-plane regions. The new colour range is 0.75 < J-K s < 2.5. 



Field 


Absolute counts (deg 2 .mag L ) 


Relative counts (%) 


1=15° b=0° 


11335 


142 


1=20° b=0° 


31734 


414 


1=26° b=0° 


18867 


131 


1=27° b=-0.5° 


24571 


150 


1=27° b=0° 


28765 


112 


1=27° b=l° 


21998 


82 



Table 3. Star count difference at the maximum of the peak in K histograms. The second column corresponds to the 



absolute counts (in number of stars per square degree per magnitude) Ndat 



N„ 



and the third one represents 



the counts (in percentage) relative to the disc Ndat 



odd 



/N„ 



jdel ■ 



3 fields have been used for the determination: 1=26° 
b=0°, 1=27° b=0° and 1=27° b=-0.5°. In these fields, the 
excess peak is rather well defined and extinction seems to 
be well modeled. The last two fields are not very reddened, 
which permits us to avoid problems of incompleteness. The 
field at 1=26° is more absorbed but the excess is clearly 
determined. Unfortunately, the fields at 1=15° and 1=20° 
are too absorbed to permit a good definition of the the 
peak excess so they cannot be used to derive the distance 
of the stars. Thus, the only fields that can be used to 
determine the distance of the excess stars are close to the 
end of the overdensity. 

The following procedure to determine the distance of 
the excess stars was used: assuming that these stars are 



of the same or similar type as the disc ones, we assigned 
as distance for each data star in a given colour-magnitude 
bin as the median of the distances of the model stars in 
the same bin. Histograms of distance distributions of the 
difference 'data minus model' are thus deduced. These his- 
tograms are presented in Figure ITTI 

The distance distribution peaks very clearly at 1=27° 
b=0°, with a maximum at 5.9 kpc. The peak is lower but 
rather well defined at 1=26° b=0° with the same max- 
imum. At 1=27° b=-0.5°, the peak is ever less well de- 
fined, with a maximum at 5.5 kpc. We can thus deduce 
from these graphs that the stars of the excesses are lo- 
cated at a distance from the Sun a little lower than 6 kpc. 
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Fig. 8. Absolute (upper) and relative (bottom) his- 
tograms of difference 'data minus model' in K s for 3 differ- 
ent hole scale length: 0.5 kpc (solid lines), 1 kpc (dashed 
lines) and 1.5 kpc (dotted lines). Units are mag^.deg -2 
for absolute histograms and % for relative ones. 
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Fig. 11. Histograms of distance distribution of the abso- 
lute difference 'data minus model'. Distance of data stars 
were derived using 12 x 12 bins of colour-magnitude with 
the following cuts: 12 < K s < 13.5 and 1.6 < J-K s < 3. for 
1=26° b=0°, and 11.7 < K s < 13. and 1.3 < J-K s < 2.5 
for 1=27° b=0°,-0.5°. 

This value confirms the result given by Hammcrslcy et al. 
(2000). 

6. Discussion 

The Besangon model has shown good performances in 
producing the observed NIR counts at intermediate lon- 
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Fig. 12. Scheme of the galactic in-plane bar 

gitudes in the Galactic plane, where the disc is the major 
contributor, even if the vertical scalelength of the model 
components may have to be fine tuned. 

We then confirm the existence of an overdensity be- 
tween 1=20° and 1=27° in the Galactic plane, which seems 
to decrease inversely to the distance from the galactic 
plane at its end (l=26°-27°), and is located at a distance 
from the Sun a little lower than 6 kpc (assuming that they 
are disc- like stars). 

The extension in longitude of the extra density and its 
confinement very close to the galactic plane may suggest 
an in-plane bar shape, as is argued in Hammersley et al. 
(2000) and Lopez-Corredoira et al. (2001). Supposing that 
it does correspond to a bar, one can estimate its bar an- 
gle and its half-length, using the distance of excess stars 
from the Sun at its top end. Neglecting the width of the 
bar end and taking a distance between the Sun and the 
galactic center R@ = 8 ± 0.5 kpc, and a distance from 
the Sun of bar stars R = 5.9 ± 0.5 kpc, we obtain (see 
Figure O that the bar has a half-length d = 3.9 ±0.4 kpc 
and an angle from the Sun-center direction 9 = 45° ± 9° . 
These results are compatible with the values obtained in 
previous works: d m 4 kpc and 6 w 43 deg (Hammersley et 
al., 2000). Moreover, other authors have given similar re- 
sults for the geometrical parameters of the Galactic bar by 
using different methods. Sevenster et al. (1999) extracted 
an angle of 44° for the bar through kinematic analysis of 
OH-IR stars. Nakai et al. (1992) studied the CO distri- 
bution in the inner Galaxy and suggested a bar near 45°, 
but with non-negligible uncertainties. More recently, Van 
Loon et al. (2003) studied luminosity distributions across 
the Galactic plane claimed for a bar with an orientation 
of 40°. 

Unfortunately, if we can estimate bar parameters as- 
suming as true the hypothesis that the excess stars be- 
long to a bar population, the fact that the study was con- 
fined at 1>15° because of the completeness of data makes 
it impossible to confirm. Even the field at 1=15°, very 
reddened and perhaps a little contaminated by the outer 
bulge, must be used with great care. Also, the most inner 
fields are dependent on the value of the disc hole scale 
length. Furthermore, the field at 1=21° b=0°, where the 
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excess disappears, shows that the structure has inhomo- 
geneities. On the other hand, extended star excesses can 
appear locally, close to rings for instance, without corre- 
sponding to a bar in term of dynamics. 

New observations are needed to settle these uncertain- 
ties. If there is a bar, an excess of stars must exist at 
negative longitudes until the other end of the hypothetic 
bar. With a half-length of 3.9 ± 0.4 kpc and an angle of 
45° ± 9°, the bar far end stars would be located at the 
longitude 1=-14° ± 2° and at a distance from the sun of 
ll.lzhl kpc. Only much deeper data will allow us to detect 
them at these longitudes, and counts will be contaminated 
by the triaxial bulge in the inner regions. Moreover, kine- 
matic measurements, at the top end of the hypothetic bar 
for instance, are necessary to know whether the extra den- 
sity corresponds to a bar in term of dynamics or not. Such 
a study, using radial velocities, is being planned. 
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